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BOUNDARY SINGULARITY OF POISSON 
AND HARMONIC BERGMAN KERNELS 


Miroslav Englis 


Abstract. We give a complete description of the boundary behaviour of the Poisson 
kernel and the harmonic Bergman kernel of a bounded domain with smooth boundary, 
which in some sense is an analogue of the similar description for the usual Bergman 
kernel on a strictly pseudoconvex domain due to Fefferman. Our main tool is the 
Boutet de Monvel calculus of pseudodifferential boundary operators, and in fact we 
describe the boundary singularity of a general potential, trace or singular Green 
operator from that calculus. 


1. Introduction 

Let Q be a bounded domain in C” and the subspace of all holomorphic 

functions in It is a standard consequence of the mean-value property of 

holomorphic functions that the point evaluations / f{z), z £ Q, are bounded 

linear functionals on and thus possesses a reproducing kernel — 

the Bergman kernel B{x,y) of Q; namely, G for each y, and 

fiy) = {f,Bi-,y))=[ f{x)B{y,x)dx V/G Lhoi(n), Vy G O. 

Jn 

If n is sufficiently nice (smoothly bounded and strictly pseudoconvex), then the cel¬ 
ebrated theorem of Fefferman [20] (with later improvements by Boutet de Monvel 
and Sjostrand [10]) gives a description of the boundary singularity of the Bergman 
kernel: namely, there exist functions a, 6 G x fl) such that 

(1) B{x, y) = + h{x, y) log p(x, y) Vx,yGlI. 

Here p{x,y) G G°°(H x H) is such that dp{x,y)/dy and dp{x,y)/dx vanish to 
infinite order on the diagonal x = y, while p{x, x) = p{x) is a defining function for 
H in the sense that p{x) > 0 for x G H and p{x) = 0, ||V/9(x)|| / 0 for x G dfl. 
(It is a consequence of the strict pseudoconvexity of H that such p(x, y) exists and 
can be chosen such that Rep(x,y) > 0 Vx,y G H, so that logp can be defined 
unambiguously.) 
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Fefferman’s expansion (1) has subsequently found far-reaching applications in 
function theory of several complex variables, complex geometry, mathematical 
physics, operator theory on function spaces, and many other areas (see [4], [21], 
[25], [27], [12] and [38] for a sample). 

The aim of this paper is to give an analogous description for the harmonic Berg¬ 
man kernel, that is, for the reproducing kernel H{x,y) of the subspace 

-^harm(^) •= {/ ^ {^) '■ f is harmouic on fl} 

on a bounded domain Q C R”, n > 2, with smooth boundary. Thus H(-,y) G 
^Lrm(^) Vy G fl and 

fiy) = {f,H{-,y))=[ f{x)H{y,x)dx yf e Ll^,^{n)yy e n. 

Jn 

(The existence of the kernel H{x,y) is again a standard, and easily verified, conse¬ 
quence of the mean-value property of harmonic functions.) 

Although harmonic Bergman kernels have been around for quite a while — 
starting probably with the book [5] — and there exist explicit formulas for them 
in a few special cases (such as the ball and the half-space [15] [2]), as well as for 
their weighted analogues [34] [30], and various applications of these in operator 
theory and function theory [14] [35] [13], the description of the boundary behaviour 
of H{x,y) in the general case seems to be lacking. To the author’s knowledge, the 
only result in this direction in the literature is due to Kang and Koo [31], who gave 
estimates for the growth of H{x,y) and its derivatives at the boundary; namely, 
for any multiindices a,/3 £ N"", 

d\°‘\+\i^\H{x,y) ^ _ 

dx^dyi^ — _l_ -g |x — y\]^+H+\l^\' 

Here d{x) = dist(x,cin) is the distance of x from the boundary. On the diagonal 
X = y, there are also analogous estimates from below. The main ingredient in their 
proof is the scaling method, familiar in the theory of several complex variables. 

The harmonic Bergman kernel is closely related to another familiar object in 
analysis, namely, to the Poisson kernel K{x,C) = K^iC), which gives the solution 
to the Dirichlet problem 

(2) A/ = 0, f\an = u 

on Q: 


(3) fix) = {u,K^)on=[ u{C)K{xX)dC, 

Jan 

where dC denotes the surface measure on dQ. The boundary behaviour of K{x,Q 
— which, again, is well known in a handful of special cases (like those mentioned 
before for H{x,y)) due to explicit formulas available — has recently been studied 
by Krantz [32], who showed (using again the scaling method) that 


K{x,C) 


Cnd{x) 
X - 


( 4 ) 
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as X and C, approach the same point on the boundary; here 


— 


r(§) 

^n/2 


is a constant depending only on the dimension n. Again, analogous estimates for 
the derivatives are also given, as well as numerous references to further information 
in the literature. 

The relationship between H[x, y) and K{x, C) is as follows. Consider the Poisson 
extension operator from (2) 


(5) K : u i-A /, A/ = 0, f\on = u, 

as an operator from L‘^{dO,) into L^(12). (In other words, K is the integral operator 
with kernel K, and K is the Schwartz (distributional) kernel of K.) It can be 
shown that K is bounded (even compact), and we denote by K* its adjoint. The 
operator 

A := K*K 

is then an injective compact operator on L‘^{dO,) with dense range, belongs to 
this range for each x G fl, and 

(6) H{x,y) = {A-^K:„Ky)an = KA-^K:,{y) Vx,yGb!. 

Phrased yet another way, H{x,y) is the distributional (Schwartz) kernel of the 
operator 

(7) KA-^K*, 


which is nothing but the orthogonal projection (harmonic Bergman projection) 
ITharm ^ L'^ {^l) OntO _ 

Our main results are the following. Let d be a smooth function on such that 
d > 0 on 0 and d(x) = dist(x, dO.) for x near the boundary, and let further x denote 
the reflection of x G 0 with respect to dQ; the latter is well defined for x near the 
boundary (see §4.3 below for the details). Finally, let 8""“^ denote the unit sphere 
in R"", R_|_ = (0, +oo) and R+ = [0, +oo). 

Theorem 1. Let 0 C R” be bounded with smooth boundary. Then there exist 
funetions F G C°°{dTl x R_|_ x 8*^“^), G G C°°{VL x 511), with = 1 for all 

( G dLl and u G 8""“^, such that 


(8) K{x,C) 


Cndjx) 
X — Ch 




+ |x - CrG'(x, C) log |x 


Cl 


for all X G Ll, ( G dLl. 

Theorem 2. Let 11 C R"" be bounded with smooth boundary. Then there exist 
functions F G C°°{Q x Q x R_|_ x 8"“^), G G G°°{Q x 0), with F{x,x,0,i') = 
n{u, Vd(x))^ — 1 for all x G did and v G 8”“^, such that 


(9) 


H{x,y) = |x-y|,-^^—f-) + G{x,y)log\x - y\ 

\x — yl"- V lx — y| / 


X - y|. 


for all x,y G il close to the boundary. 
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Examples are also given showing that, as with the original Fefferman expansion 
(1) in the holomorphic case, the logarithmic term in (8) and (9) is in general nonzero, 
even though it is absent in some special situations (e.g. for n = 2 in Theorem 1). 

We remark that the statement of Theorem 1 is equivalent to the same but with 
F(x, |x — Cl, ), F G X R_|_ X depending on x instead of C, or even 

with F(x,C, \x — Cl, ), F G C°°{0, x dQ x R_|_ x 8"'“^), allowed to depend on 
both variables (with E(C,C,0, = 1 for all C £ dQ and u G 8”“^). Similarly, one 

can replace the F{x,y, \x — y\, ZZ^-\ ) Theorem 2 by either F{x, \x — y\, ZZ^-\ ) or 

Note that the leading order terms in (8) and (9) recover, of course, the coarser 
estimates from [32] and [31], respectively, recalled above. 

As already mentioned, the main tool for deriving the estimates in [32] and [31] 
was the scaling method; in [32] there is also a rough sketch of a possible way 
how to obtain (4) — and, in the present author’s opinion, also (8) and (9) — 
by means of a microlocal reduction to a model case using Fourier integral operators 
(see e.g. §8 in [9], or [10], for a sample of this kind of techniques). We employ 
a different approach here (avoiding, in particular, the sophisticated machinery of 
Fourier integral operators), relying on the extension, originating with the work 
of Boutet de Monvel and others, of the standard calculus of pseudodifferential 
operators to boundary value problems (“boundary TDOs”). This extension includes 
the “potential” operators like our K, as well as “trace” operators like K*, TDOs 
on dQ (which turns out to be the case of the operator A = K*K above), and, 
finally, the so-called “singular Green operators”, of which an example is precisely 
the harmonic Bergman projection Ilharm = KA“^K*. The question of identifying 
the singularities of the distributional kernels of all these operators (including thus, 
among others, the distributional kernels K{x,Q of K and H{x,y) of KA“^K*) 
can then be handled — as in the case of ordinary TDOs, where such results go back 
to Calderon [11] and others — with the aid of the well-known standard apparatus 
of homogeneous distributions. In particular, in addition to (8) and (9), we in 
fact obtain a description of the singularities of the distributional kernels of general 
potential, trace and singular Green operators from Boutet de Monvel’s calculus. 
This includes, by the way, also weighted harmonic Bergman kernels R^(x,?/) on 
n with respect to weights w G (7°“ (12) that have “the same order of vanishing” 
at all points of the boundary (i.e. w{x) = d{x)^g{x), where d is as above, m > — 1 
and g G (7°° (12) is positive on 12), as well as reproducing kernels of some harmonic 
Sobolev spaces; see Section 7 below for the details. 

It should be noted that all the methods used here are well-established and, 
in fact, the real-analytic case is more or less treated in [7] (with hints about the 
C°° situation in §7 there); however, the final results seem not to be available 
anywhere in the literature, nor to be a part of common knowledge, which was the 
author’s reason for writing this paper. Since the paper hopes to aim also at some 
audience from operator theory and complex function theory, we also decided to 
include a fairly more extensive review of the various prerequisites than would be 
strictly necessary for experts in the area of TDOs. 

The paper is organized as follows. Section 2 reviews the necessary material on 
boundary 'kDOs (Boutet de Monvel’s calculus). Various facts about the Poisson 
kernel and the harmonic Bergman kernel are collected in Section 3. Singularity 
resolution of the Schwartz kernel of a general boundary TDO is given in Section 4. 
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Some more specific calculations, which allow us to refine these results to obtain 
Theorem 1 and Theorem 2, are presented in Section 5 and Section 6, respectively. 
The final Section 7 contains miscellaneous concluding remarks, open problems, and 
the like. 

Notation. Throughout the paper, Z and N stand for the integers and the 
nonnegative integers, respectively. For a multiindex a G N"", 5“ (or just 5“ if 
the variable is clear from the context) stands for the partial derivative /dx°' = 
on R"", and = (—The Euclidean norm of x G R"' 
is denoted by ||x|| or simply even |x|; the Euclidean scalar product of x,y G R”^ is 
denoted by (x,y) or x • y. The inner products in L^(n) and L‘^{dVL) are denoted 
by (•, •)o and (•, •)ao, respectively, with the subscript omitted if there is no danger 
of confusion. Finally, unless explicitly stated otherwise, fl C R"" will be a bounded 
domain with smooth (=(7“) boundary. 

2. BOUTET DE MONVEL CALCULUS 

2.1 PseudodifFerential operators. Recall that a pseudodifferential operator 
(TDO for short) on R"" is an operator of the form 

(10) Au{x) = f e^^'^a{x,^) 

JR" 

Here = (27r)“"'d^ is a renormalization of the Lebesgue measure d^, x ■ ^ = 
H-^ Xnin, and 


(11) II(0 = / u{x)e dx 

JR" 

is the Fourier transform of u. The operator (10) is also written as a{x,D), a no¬ 
tation justified by the fact that A = f{x)D°‘ for a(x,^) = /(x)^“. The function 
a(x,^), called the total symbol of A, is usually assumed to lie in Hormander’s class 
5™'(R”,R”) for some m G R, i.e. to belong to C'°°(R"' x R”) and satisfy 


( 12 ) 


19) 


'dfa(x,OI 


a:eA,^6R" (1 + ICI)"" 


< OO 


for any compact subset K of R"" and any multiindices a, (3. It is standard that 
the operator a(x, D) is then well-defined for any u in, say, the Schwartz space 
5(R”), and extends by continuity to a bounded operator from the Sobolev space 
lF<f„^p(R-) into K;”*(R"), for any s G R, as well as from V{K^) into C'°°(R’") 
and from the space P(.Qj^p(R"') of compactly supported distributions into V'CRA). 
Most TDOs in this paper will be classical (or polyhomogeneous), meaning that 
their total symbol has an asymptotic expansion 

OO 

(13) a{x,^) 

1=0 

where am-j is C°° in x,^ and positive homogeneous of degree m — j in for 
|.^| > 1. The notation means, by definition, that the difference between a and 
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X]J=o should belong to the Hormander class S^~^, for each k = 0,1,2,.... 

The set of all classical TDOs as above (i.e. of order m) will be denoted by the 
(larger) class of all (not necessarily classical) TDOs with total symbol in S'™' will 
be denoted by T™, and we set, as usual, Td := UmeR'^d’ 

:= PlmeR ~ rimeR classes T™ make sense even for complex m, 

but we will not need those in this paper.) The operators in are smoothing 

operators, i.e. their Schwartz kernel is in C'°“(R"' x R"). 

Combining (10) and (11), we can write the definition of a{x,D) = ^ as 

(14) Au{x) = jj e^^^~y^'^a{x,^)u{y)dyd^, 

with the double integral interpreted in suitable sense. One can then define a more 
general variant of TDOs A = a{x, y, D) allowing the symbol a (called then an 
^^amplitude”) to depend also on y: 

Au(x) = jj e"^'^~^^'^a{x,y,()u{y)dyd^, 

where a G C°°(R^’^) is again assumed to satisfy (12) for some m G R (with 5“ 
replaced by d^d'^, and the supremum extended also over all y G iiT). It turns out, 
nonetheless, that one obtains essentially the same class of operators in this way: 
namely, for any amplitude a{x,y,^) G S'™(R^”',R*^) there exists a symbol b{x,^) 
in S'™(R”^,R”^) such that a{x,y,D) — b{x,D) is a smoothing operator. Modulo 
smoothing operators, 6(x, y) is uniquely determined by the asymptotic expansion 


where the summation extends over all multiindices a. 

A 'kDO A = a{x, y, D) is properly supported if both A and A* = a{y, x, D) 
(the bar stands for complex conjugation) have the following property: for each 
compact set K C R"", there exists a compact K' C R”' such that distributions 
supported in K are mapped into distributions supported in K'. Every TDO can 
be written as the sum of a properly supported TDO and a smoothing operator. 

In general, TDOs do not preserve the support of a function (or distribution), 
however, they have the pseudolocality property: if [/ C R"" is an open set and u is 
C°° on U, then Au is also C°° on U. 

The composition of two (classical) TDOs, at least one of which is properly sup¬ 
ported, is again a (classical) TDO: if a G T™, b G we have a{x,D)b{x,D) = 
c{x,D) modulo smoothing operators, where c G satisfies 

(15) c(x, 0 ~ y] —.D1a{x, 05“6(x, 0- 


(The right-hand side determines c only modulo symbols in S~°°.) 

TDOs behave nicely under coordinate transformations and, consequently, one can 
define the classes T™, T™ also on manifolds. (Namely, for any partitions of unity 
{i^j} subordinate to some atlas {4*j} of local coordinate charts, with ipj = 1 
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on supp(/>j, one declares an operator A on the manifold to belong to etc. if and 
only if the operators 

(16) 

are in etc. for all j,k. Note that the operators (16) are automatically 

properly supported.) In particular, we have the classes 'l')y(clfl), m G R, of clas¬ 
sical 'I'DOs on the boundary dQ of a smoothly bounded domain Q. Note that 
an operator A G 'h)j(9fl) is then continuous from the Sobolev space iy<foj„p(cin) 
into IV[™“®(cin), for any s G R; in particular, if Q is bounded, A G 'I')l|^(9n) is 
continuous from hF^(9n) into 

Throughout the rest of this paper, we will almost exclusively be concerned with 
classical TDOs on flfl, and we will thus usually abbreviate lT®(cin), T"*(fin), etc., 
just to IT®, 'h™', and so on, and omit the adjective “classical”. 

Standard references for the material above are e.g. the books by Hormander [29], 
Treves [41], Shubin [40], Tolland [23] or Grubb [26]. (All these give also a treatment 
of TDOs on the domain Q C R”, rather than on its boundary dQ; this however 
will not be needed here.) 

2.2 Boundary TDOs. There is an extension of the theory from §2.1 to manifolds 
with boundary, which we now review; a good reference for the material below and 
in §§2.3 and 2.4 is Rempel and Schulze [37] and Grubb [26]. 

Gonsider first the case of the upper half-space R” = R""”^ x R_|_ in R”, R!J: = 
{x = {x',Xn) G R"""^ X R : > 0}. Let us denote by 5+ the space of restrictions 

to R+ = [0,-|-oo) of functions in the Schwartz space 5(R), and, similarly, by 
the space of restrictions to R_|_ x R_|_ of functions in the Schwartz space 5(R^). For 
d G R, the space 5^(R’”, 5+) consists, by definition, of all C°° functions / 

of {X,Xn,i') £ R-™" R-i- lying in 5+ with respect to Xn, such that for all 

a G /3 G N™", A:, Z G N and compact subsets K of R™", 

\XndLdx9pf{X,Xn,^')\ 

sup —;- , T , , , —rn— < oo. 

xeA,.„>o (1 + 

The space 5^(R"^, R^-^ 5++) consists similarly, by definition, of all C°° functions 
g of {X, Xn,yn,C') £ R-"* X R-i- X R_(_ X lying in with respect to {xn,yn), 

such that for all a G /3 G N™", /c, Z,p, g G N and compact subsets K of R™", 

\xldiy^dU^x^<^,g{X, Xn,2/n,?0l 

The subspaces of classical (or polyhomogeneous) elements in these S'^ are defined 
to consist, analogously as for symbols of TDOs in §2.1, of all f,g with asymptotic 
expansions 

OO OO 

or gr^'^gd-u 
1=0 z=o 

respectively, where 

f,{X,^,Xe)=X + ^fiX,Xr.,a, gj{X,^,l^,Xa = X^^^9iX,Xn,yn,e), 
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for A > 1 and |.^'| > 1. 

A negligible potential operator is an operator from functions on into func¬ 

tions on R" whose distributional (Schwartz) kernel is smooth on R” x 
A potential operator of order d G R is an operator from functions on R”^“^ into 
functions on R!J: which is a sum of a negligible potential operator and an operator 
defined by 

(17) Ku{x',Xn)= f k{x',Xn,i')u{i')di' 

where k G S'^[~^(R”^“^,R"'“^,5_|-). As with 'I'DO, one gets the same class of opera¬ 
tors using kernels k{y',Xn,C') instead of k{x',Xn,C')j and even symbols (amplitudes) 
depending on both x' and y' can be allowed: 

Ku{x',Xn)= [ [ k{x',y',Xn,^')u{y')dy'd^', 

with k G 5;ff^(R2^-2,R^-i,5+) and the integral suitably interpreted (as an oscil¬ 
latory integral), but again yield the same class of operators. 

A negligible trace operator of class r G N is an operator from functions on R!J: 
into functions on R""”^ of the form 

Sjlj + T, 

0<j<r 


where 'jj is the j-th order boundary normal derivative ^jU^x') = (d^^u)(x',0), and 
Sj, T have distributional kernels that are smooth on R”“^ x R"“^ and R*^"^ x RIJ:, 
respectively. A trace operator of class r G N and order d G R is an operator from 
functions on R!J: into functions on R”“^ which is a sum of a negligible trace operator 
and an operator of the form 

(18) ^ Sj^j + T, 

0<j<r 

where Sj G (R"""^) and 

p poo 

Tu{x')= / / t{x',Xn,i')u{i',Xn)dXndi', 

dRr*-i do 

where u{i'■, Xn) stands for the partial Fourier transform of u{x', x„) with respect to 
the x' variable, and t G 5^j(R”“^, R"""^, 5+). Again, symbols t{y' ,Xn,0 or more 
general amplitudes t{x',y',Xn,C') from 5)?[(R^’^“^, 5+) can also be allowed, 

but lead to the same class of operators. 

A negligible singular Green operator of class r G N is an operator on R” of the 
form 

+ G, 

0<j<r 

where Kj are negligible potential operators and the distributional kernel of G is 
smooth on R” x Rlj:. A singular Green operator of class r G N and order d G R is 
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an operator on R” which is a sum of a negligible singular Green operator and an 
operator of the form 

(19) K,^j+G, 

0<j<r 

where Kj are potential operators of order d — j and 

(20) Gu{x',Xn)=[ [ g{x',Xn,yn,i')u{i\yn)dyndi', 

90 

whereof G 5++). Again, symbols g{y'■, Xn,yn, CO more general 

amplitudes g{x',y',Xn,yn,^') from S'^j“^(R”'“^ x R"'“^,R”“^,5_(__(_) can be allowed 
but lead to the same class of operators. 

We denote the three types of operators just introduced by /C'^(R!j:), 70*(R!(), 
and 0(?(R!j:), respectively. 

As with 'hDOs, the three types of operators above are pseudolocal in the x'- 
variable, in the obvious sense. They are called ^^boundary TDOs” (bTDOs). 

We will call the functions k,t,g the symbols of K, T and G, respectively.^ 

The calculus of bTDOs, initiated by Boutet de Monvel [6], [8], implies, first 
of all, that the three classes of operators defined above behave nicely under compo¬ 
sitions and taking adjoints; and, second, that solution operators to boundary value 
problems — like our Poisson operator K which we had before — belong to this 
calculus. Let us briefly describe only those details about this that we will need. 

Remark. There exist also more general “nonclassical” versions of T, K and G, but 
they will not be needed in this paper. □ 

2.3 Composition rules. Let iL be a potential operator, T a trace operator of 
class zero, and S a (classical) TDO on R""”^. Then 

KS is a potential operator, 

ST is a trace operator, 

TK is a TDO on dVl, and 
KT is a singular Green operator, 

with orders adding up, provided at least one of the factors in each case is properly 
supported in the x'-variable. The symbols of these products are given by the usual 
formula (15) in the x^variable, but new things happen in the Xn-variable. Let us 
define 

{k t){x', Xn,yn,^') ■= k{x', Xn,Ot{x\ 2/n, ?')> 

t{x' ,Xn,C')k{x' ,Xn,^') dXn, 

{kOn s){x',Xn,C') ■= k{x',Xn,(')s{x',C'), 

{sOn t){x',Xn,C') ■= s{x',C')t{x',Xr„C'). 

^The reader is warned that this is somewhat at odds with the standard notations and termi¬ 
nology in the area, including the references such as [26] mentioned further down: namely, our k is 
usually denoted by k and called (for good reasons) symbol-kernel rather than symbol, of K € 
the term “symbol” being reserved for the Fourier transform of our k{x’ with respect to 
the Xn variable, after extending by zero for Xn < 0. Similarly for T and Q. Also “potential 
operators” are called “Poisson operators” in [26] and elsewhere; we reserve here the latter term 
for the Poisson operator K from (5). 


( 22 ) 


{tOn k){x',i') := [ 
Jo 
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Let {A, B) be one of the pairs of operators {K, S), {S, T), (T, K) or (iL, T), and let 
a, b be the corresponding symbols k, t, s or g. Then AB = C is an operator of the 
type indicated in ( 21 ), with symbol 

(23) c~^^(Zl5“,a)o„(9:,6). 

q; 

There are also composition formulas for trace operators of class r > 0; we will only 
need the case of TK with r = 1 here, which reduces to the simple rule 

(24) = S with s{x' ,^') = k{x' 

Similarly for 70 G with a singular Green operator G. Finally, there are rules for 
the products PK = K' and PG = G\ where P is a differential operator on R"": 
namely, on the level of the x' variables one uses the standard TDO rule (15), while 
for P = x^d^.^ one just has 

k'{x',xn,i') = x^dl^k{x',xn,^') 
and similarly for g and g'. 

Similar rules hold also for other possible compositions, but again will not be 
needed here. 

In addition to composition, the calculus behaves well with respect to adjoints: 
namely, the adjoint of a Poisson operator K G JC‘^ with symbol k is the trace 
operator T G with symbol 

(25) t{x',xn,i') ~ V ^T»7<9“,A:(P,Xn,C') 

ai ^ 


(the bar denoting complex conjugation). There is also an analogue for Qq, which 
will not be needed here. 

2.4 Boundary value problems. For T G 77”(R+), the boundary symbol of T 
is the operator t{x',^',Dn) from 5_|_ into C defined by 

pCO 

t{x',C,Dn)u= t{x',Xn,C)'uiXn)dXn+ ^ Sj{x',C)7jU. 

0<j<r 

Similarly, the boundary symbol of G G is the operator on 5_|_ defined by 

pOO 

g{x',^',Dn)u{Xn)= g{x',Xn,yn,^')u{yn)dyn + ^ kj{x',Xn,^'hjU, 

0<j<r 

and the boundary symbol of iL G (R” ) is the operator from C into 5-(_ defined by 
k{x', Dn)a = k{x', x^, ^')a, a G C. 

The principal boundary symbol is defined similarly using only the leading-order 
term (the principal symbol) tm, Pd or kq in the asymptotic expansion of t, g or k, 
respectively. 
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Now let A be the matrix of operators 


(26) 


A = 


F+ + G 
T 


K 

S ’ 


with K, T, G as above, S a 'I'DO on of order m+q—d, and = r_|_Pe_|_ where 
P is a differential operator on R"", of the same order d as G, with smooth coefficients, 
r+ stands for the restriction from R” to R!(, and e+ for the operator of extension 
by zero from R” to R”; abusing notation slightly, we will sometimes write just P 
instead of Also, more generally we even allow K to be an M-tuple (row matrix) 
of potential operators, T an M'-tuple (column matrix) of trace operators, and S 
an M' X M matrix of ^DOs (M, M' = 0,1, 2,...). The principal interior symbol 
of A is just the principal symbol Pd{x,^) of P; the principal boundary symbol is 
the matrix of operators 




Pd{x',0,C',Dn) + gd{x', , D^) 

tmix',C',Dn) 


kq{x',C',Dn) 
Sm+q—d{x , ^ ) 


with Pd{x', 0 ,^', Dn) ; 5+ —)■ 5+ defined by 


Pd{x',0,^',Dn)u{Xn) = / e*^'*^"Pd(x',0,^',^n)^(^n)(?^n, 

Jo 

where e+ is the operator of “extension by zero” from R+ to R. The matrix A is 
said to be elliptic (of orders d,m,q and class r) if Pd{x,^) / 0 for all |^| > 1 and 
all X G R”, and a{x',^',Dn) ■ 5+ x — )■ 5+ x is bijective for all |^'| > 1 

and x' G R”“^. 

Finally, the three classes of operators JC,P,Q are invariant under coordinate 
transformations of R” preserving the boundary 9R!J: = R""”^ x {0}. Thanks to 
this invariance, one can again — via local coordinate charts — dehne the analogous 
classes /C'^, and for R” replaced by any manifold with boundary; in particu¬ 
lar, the classes }C^{Q), T^{^) and Gf{^) for R” replaced by a domain Q. C R”, and 
the boundary cfR” = R*^"^ replaced by the boundary dGl of Gl. The dehnitions of 
principal boundary symbols and ellipticity carry over to this setting. 

With all these definitions, it is then true that 


if A is elliptic of orders d, m, q and class r and Q, is bounded, 
then there exists a matrix B of operators of the form (26), ellip¬ 
tic of orders —d,—q,—m and class r' = max{r — d, 0 }, which is 
a parametrix for A, i.e. AB — I and BA — I are negligible. 


A potential operator K G /C'^(ff) whose amplitude is compactly supported in 
both x' and y' can be shown to map VF®(5n) continuously into (ff), for 

any s G R. In particular, if Q. is bounded, then any K G /C^ is continuous from 
VF®(5n) into s G R. Similarly, T G T^{Gl) and G G Qf{Gi) are 

continuous IF'*(n) —)■ IT®“‘^“^/^(cA 2 ) and IT®(ff) —)■ IT^“'^(n), respectively, for 
any s > r — i, s > 0 , if If is bounded. 

Remark. The result (27) actually holds also for suitable pseudodifferential (i.e. not 
necessarily differential) operators P, namely, for those that satisfy the transmission 
condition; but again we will not need this here. □ 
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Remark. For unbounded, (27) gets more complicated; for instance, even for 
the matrix (26) corresponding to the ordinary Dirichlet problem for A on R” 
(cf. the beginning of Section 3), the kernel k{x',Xn,C') in (17) of the corresponding 
Poisson operator K equals k{x',Xn,C') = I, which is not smooth at = 0. 

Nonetheless, the localized operator Kx (i.e. K preceded by multiplication by x); 
with any x £ 2?(R"“^), will already be an honest potential operator, i.e. with 
smooth k. □ 


3. Some simple facts 

In this section, we establish various useful facts about the Poisson operator K, 
the Poisson kernel K{x,C), the harmonic Bergman projection Ilharm) and the har¬ 
monic Bergman kernel H{x,y), that will be needed later. 

Consider first our Poisson (harmonic extension) operator K : n i—)■ / solving the 
classical harmonic extension problem 

A/ = 0 on fl, f = u on dk}; 


that is. 


The operator 


A 

7o 


/ 






C“(n) 

© 

C°°{dQ) 


is of the type (26) discussed in §2.4, with interior symbol P 2 {x,^) 
boundary symbol given by 


a{x',^',Dn) = 


7o 




5 + 

© 

C 


— and 


has the unique solution u{xn) = ae I in 

is uniquely solvable. Thus A is elliptic and has 

a parametrix, unique up to negligible operators, of the form [P + G,K], with P 
a TDO on R” of order — 2, G a singular Green operator of order —2 and class 
max{l — 2,0} = 0, and K a potential operator on n of order 0. On the other hand, 
we know that 

R = [G,K], 

with K as above and G the solution operator to the classical Dirichlet problem 


The equation a(x', Dn)u = 
and similarly a{x',Dn)u = 


AG/ = / on D, Gf = 0 on dkl, 


is a right inverse for A. It follows that we can take P + G = G and K = iF; that is, 
(28) KG/C°(n). 
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By the facts reviewed in Section 2, it follows that 

(29) K:W^{dn) Vs G R; 

in particular, K is a continuous operator from W^{dVL) = L‘^{dVL) into C 

L^arm(^) {^) (even compact). Let K* : L‘^{VL) —)■ L‘^{dVl) be its Hilbert space 
adjoint and consider the operator 


A := K*K. 

Note that A is injective (since K is). From 

(30) A“^K*K = I 

it is immediate that KA“^K* is selfadjoint, idempotent, and vanishes on KerK* = 
(RanK)-*- while being the identity on the range of K. Consequently, 

(31) KA-lR* = Hharm 

is the orthogonal projection in L^(n) onto the closure of RanK, i.e. the harmonic 
Bergman projection of L^(fl) onto Lharm(^)- 

Remark. Comparing (30) with the definition of K, we also see that 

A“^K*|RanK =: 7 

coincides with the operator of “taking the nontangential boundary values” of har¬ 
monic functions. By elliptic regularity (or Boutet de Monvel’s calculus), 7 again 
extends to a continuous operator from IFharm(^) for any s G R, 

which is the left inverse of K. □ 

By Boutet de Monvel’s calculus, we have from (28) K* G 7j)~^(f2) and thus 

A = K*K G 

Furthermore, by ( 22 ), (23) and (25), the leading symbol s_i of A is given by, in any 
local coordinate chart. 


S_l(x',^') = / \ko{x\Xn,i')\^ dXn> 

Jo 

i.e. A is elliptic. Rs inverse A“^ is thus a positive injective selfadjoint operator in 
'I'Jj(cA2). Finally, by Boutet de Monvel’s calculus again and (31), 

Hharm = KA'^K* G 


is a singular Green operator on kl of class zero and order 0. 

We conclude this section by proving one more formula relating the Poisson and 
the harmonic Bergman kernel. 
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Proposition 3. For x,y £ Q, 

H{x,y) = {A-^K^,Ky)9n. 


Proof. For any / G Lf{Pt) and x we have 


Il-harm/ {x) 



f{y)H{x,y)dy. 


Thus 


= nharm/(x) = KA-^K* f {x) 
= {A-^K*f,K^)9n by (3) 
= {f,KA-^K:,)n. 

This means that and 

H^{y) = ^A-^K^{y) = {A-^K^,Ky)9n 


as claimed. □ 


4. Schwartz kernels of boundary TDOs 

4.1 Homogeneous distributions. Recall that a function u G C'°°(R"' \ {0}) is 
said to be homogeneous of degree s, s G R, if u{Xx) = A®u(x) for all x G R” \ {0} 
and A > 0. The definition extends also to distributions on R"^ \ {0} or R” in an 
obvious way. It is standard (see e.g. [28], Chapter 3, §2) that u can always be 
prolonged from R"" \{0} to a tempered distribution u on R"" and, furthermore, ii is 
uniquely determined and homogeneous if — s — n ^ N, while for s = —n — k, k G N, 
ii is not unique (one can add any derivative of order k of the Dirac distribution 6 
at the origin) and is not homogeneous in general, but contains a logarithmic term. 
The Fourier transform of ii has the form 

(32) ii = U - Q\og \ ■ \ , 

where U, U £ C'°°(R"' \ {0}) if u G C'°“(R’^ \ {0}), is a homogeneous distribution 
on R"" of degree — n — s, while Q is a polynomial of degree — n — s =: k (thus <5 = 0 
if /c ^ N) given by 

( 33 ) Q(i)=f thipAUd., 

where dx stands for the unnormalized surface measure on the unit sphere 

These facts can be neatly summarized as follows: for s G R, let FLs be the vector 
space of all distributions on R” of the form 

Us + Qs log I • I + P-n-siD)6o, 

where Us £ C'°“(R"'\{0}) is homogeneous of degree s, Qs is a polynomial of degree s 
(thus Qs = 0 if s is not a nonnegative integer), P-n-s is a polynomial of degree 
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—n — s (thus P-n-s = 0 if —n — s is not a nonnegative integer), and stands 
for the Dirac distribution (unit point mass) at the origin. Then the prolongation ii 
of a function u G C'°“(R" \ {0}) homogeneous of degree s always exists, is unique 
if —n — s ^ N, and all such prolongations belong to T-Lg- Noting that the Fourier 
transform of log |x| coincides (up to a constant factor) with it further follows 

that in fact 

(34) the Fourier transform is a bijection of Jia onto P-n-s, Vs G R. 

The following assertions seem to be standard knowledge, but difficult to pin¬ 
point in the literature in complete generality; see e.g. Neri [36], Theorem 1.5, 
or Calderon [11], Theorem 28. 

Recall that a function y G T>(R"), 0 < X < 1) is called a cutoff function if x = 1 
in some neighbourhood of the origin; and 0 is a patch function ii 6 = 1 — X; with x 
a cutoff function. 

Proposition 4. (i) Let u G (^“(R” \ {0}) he homogeneous of degree s, s G R, 
and let 6 be a patch function. Then 

(35) eu= {U - Q\og\-\)x +g, 

where U and Q are as in (32) and (33), x is a cutoff function, and g G 5(R”). 

(ii) Let 9, x be as above. For s G R, denote 

9Ps + S : = {9u + g : u £ Ps,g G 5}, 
xPs +S : = {xu + g : u GPa,g £ S} 

{these sets are independent of the choice of x o,nd 6). Then the Fourier transform 
is a bijection of OPg + S onto xP-n-s + S. 

Proof, (i) Let u be as above; then 6u = 9u differs from h by a distribution with 
compact support. By the Paley-Wiener theorem, we thus have 

9u = Oil = ii + h, 

with h an entire function on C”. Hence, for any cutoff function X; Gu — xu £ 
C'°°(R"). 

On the other hand, for any multiindices a,/?, = (—l)l^lD“x^0tt is the 

Fourier transform of a smooth function on R"" homogeneous of degree s + |/3| — [al 
near infinity, and, hence, integrable for |a| > s + \/3\ + n. Thus by the Riemann- 
Lebesgue lemma, ^°‘Di^9u £ C'o(R"') for |q;| > s + |/1| -|- n. It follows that d^^Ou — 
and, hence, d^{9u — xu) — decays rapidly at infinity, for any /3. Thus 9u — xu £ 
5(R”), proving the first claim. 

(ii) Just as in the proof of part (i), one sees that 6u — xu £ S ior u £ Pa (the only 
difference is that D°'x^9u need now no longer be homogeneous of degree s-|- 1 / 1 | — |a| 
near infinity, but may also contain a term of this form multiplied by log |x|, which 
however has no effect on its integrability). Thus by (32), 6Pa C xP-n-s + 5, and 
hence also [9Pa + 5]^C xP-n-s + S. 

Conversely, if u G P-n-a, then by (34) there is m G Pa with v = u, so as we 
have observed at the beginning of the preceding paragraph, 6u — xv £ 5; thus 
Xi; G + 5. Hence xP-n-s + S C ^a + S = [OPa + 5]T □ 
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Corollary 5. If A G is properly supported, then its Schwartz kernel Ia 

satisfies for x ^ y 

kA{x,y) = |x - |x - y|, + G(x, y) log |x - y|, 

where F G C'“(R’" x R^ x G G C'°“(R2’"), ifmGZ,m> -n; 

kA{x,y) = G{x,y) log |x - y| + f(^x,\x -y\, 

where F G C'°°(R"' x R_|_ x 8 ""“^) and G G C'°°(R^”') vanishes to order —n — m at 
X = y, if m G Z, m < —n; and 

kA{x,y) = |x-yr’"“”"F(^x, lx-y|,||—+G(x,y), 
where F G C'°°(R"' x R+ x 8 "'“^), G G C'°“(R^"'), i/m ^ Z. 

Proof. This follows from the last proposition by a kind of argument which is very 
standard; we put it down in some detail here for the record but will be more brief 
on similar occasions later on. 

In view of (14), the Schwartz kernel of A is given by kA{x,y) = a(x,x — y), 
where d{x,z) = F^ 2 a(x,^) is the inverse Fourier transform of the symbol a oi A 
with respect to the / variable. Using the polyhomogeneous expansion (13) of a, 
we can write, for any N = 0,1, 2,... , 


x-y \ 

X - y|/’ 


N-l 

a{x,0 = ^(0 ^ am-j{x,0 + a(N){x,0, 

i=o 

where 0 is a patch function, Om-j G C'“(R”' x (R"" \ { 0 })) is homogeneous in of 
degree m — j, and a(N) G . Using the preceding proposition, this yields 


N-l 

dix,z) = xiz) E E — m — n (x, z) + Qj-m-nix, z) log jz 
J=0 


+ g(x,z) + a(Ar)(x,z), 


with y, Uj-m-n, Qj-m-n, 9 as in the proposition, the latter three depending 
smoothly on x. By homogeneity, we can write 


Gj — jri — 7 


{x,z) = \z\^ ”Uj_m_n(x, |f|) 


with Uj-m-n £ C'°°(R’^ X 8 "“^). A standard argument (imitating the proof of 
Borel’s theorem, see e.g. [40], Proposition 3.5, or [29], Proposition 18.1.3) produces 
a function il(x,r, C), il G C'°°(R"' x R X 8 "' ^), whose /-th order term in the Taylor 
expansion in the second variable is Fu{x,f): that is. 


N-l 


il(x, r, C) = ^ Fu(x, C) + 0{r^) 
i=o 
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for any = 0,1, 2,... . Similarly, one can find O G such that 

m+n+W —1 

0(x,2;)= ^ Qj-„^-n{x,z) + 0{\z\^) 

j=max{ 0 , m+n} 

for any N = 0,1,2,.... Since a(Ar) G UN > m + n + 1 

(see e.g. [23], Theorem 8.8), and N = 0,1, 2,... can be taken arbitrary, it follows 
that the function 

c{x, z) = a{x, z) — \z\, ■^) — £l(x, z) log \z\ 

belongs to £qj, hence, to C'°°(R^”'). Now £2 = 0 if m ^ Z, 

while for m G Z the function c{x, z) can be combined together with the middle 
term on the right into \z\~'^~'^F{x, \z\, j^), and the assertion follows. □ 

Note that since x—y is a smooth function of \x—y\ and ]f5^) we get an equivalent 
statement if we replace F{x,\x—y\, ) in the last corollary hy F{y,\x—y\, , 

or even F{x, y, \x — y\, |f5||) with F G C'“(R^"' x R_|_ x 

The last proof in fact establishes (omitting the variable x, which features only 
as a smooth parameter throughout) the validity of the first part of the following 
proposition; the second part is proved in the same way. Let z = r( {r G R+,C £ 
Sn-i) i-jg polar coordinates in R"". For s G R, let us denote by Cs the vector 
space of all distributions on R"" of the form 

(r*F(r,C))- +G(rC)logr + P(T>),5o + R, 

where F,G vanish for r > 1, F G C'°“(R+ x S"“^), H G 5(R"'), P is a polynomial 
of degree —s — n (thus P = 0 if —s — n ^ N), and G G C'°“(R”) vanishes to order s 
at the origin if s G N, while G = 0 is s ^ N. Here the prolongation (r'^P)’ of r®P 
from R” \ {0} to R” exists, since one can write the latter, choosing A: G N so large 
that s + A: > —n, as the sum 


PP(r,C) = 


fc-i 

E' 

j=0 


s+j 


5tgL)+r-+‘F#(r.C) 


of homogeneous functions of degrees s,s+l,...,s + A: — 1 and a function of moderate 
growth in Lj'Q^(R”). Finally, let Vm denote the space of all polynomials of degree 
< m on R"", and let again x be a cutoff function and 6 a patch function. 

Proposition 6. (i) The Fourier transform T is maps the Hormander elass S'^, 
m G R, into C-n-m- 

(ii) In faet, F is a bijeetion of 


+ eVm log := {a + 6plog \ ■ \ : a G G Vm} 


onto C—Yi — m' 

Proof, (i) By definition, 5™" C 5 + (in the obvious sense, i.e. for 

/ G 5™" there exist hj G TLm-j such that / — hj G for all G N; 
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in fact, the hj furthermore do not contain any logarithmic terms, nor the derivatives 
of 5o which are killed by 9). Using Proposition 4(i) and the same argument as in 
the last proof, this implies C 5 + = C-n-m- 

(ii) Now 5*” + OVm^og = 5 + (this time we use also the log¬ 

arithmic terms in Tim-j, while the derivatives of ho still play no role); since 
C-n-m = 5 -|- xYlij>o'^-n-m+j-, the claim follows in the same way using part 
(ii) of Proposition 4. □ 

Using local charts, one gets also the analogue of Corollary 5 for operators on dCl. 

Corollary 7. // 14 C R” is a bounded domain with smooth boundary and A G 
'I'^((914), then its Sehwartz kernel {with respect to some smooth volume element 
on dQ) is of the same form as in Corollary 5, with F G C°°{dCl x R+ x S”“^). 

Proof. Immediate from Corollary 5 and the observation that, for any diffeomor- 
phism 4>, and are smooth functions of x, \x-y\ and |f 5 f| • □ 

Let us now give the same treatment to boundary 'kDOs. 

4.2 Potential operators. We start with the Poisson type operators. 

Theorem 8. If K £ JC‘^{Il), where 14 C R"" is a bounded domain with smooth 
boundary, or Cl = R”, then its Schwartz kernel kx satisfies 

kxix, C) = |x - |x - Cl, |^_^| ) + G{x, C) log |x - Cl 

if d £ 7i, d > 1 — n; 

kxix, C) = G{x, C) log |x - Cl + i^(C, k - Cl, 
with G{x, C) vanishing to order 1 — n — donx = C,ifd£Z, d<l — n; and 
kxix, C) = |x - Cl^“’"“'^-?^(C, k - Cl, ||^ _^| ) + G{x, C) 

i/d ^ Z, /or X G 14, C e dCl, where F £ C°°{dCl x ^ x G £ G°°{Ti x 514). 

Note that by Seeley’s extension theorem, we can even assume that F £ C'°“(R”^ x 
R X R*^) and G £ C'°°(R" x R*"). 

Proof. It is enough to deal with the case of 14 = R” and K of the form (17) 
(the negligible piece of the operator contributes a Schwartz kernel smooth on R!j: x 
R!J:), where it is now convenient to take the y'-version of the kernel k, i.e. 

Ku{x) = [ [ e'G'-y’X’k{y',Xn,f,')u{y')dy'df,' 

interpreted as oscillatory integral, with k £ 5)?j~^(R"'“^,R”^“^,5_(_). The Schwartz 
kernel is thus given by 


kxix, y') = k{y', Xn,x' - y') 
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where k denotes the inverse Fourier transform of k{y', Xn,^') with respect to 
We claim, first of all, that 

(36) k G ^ ke X if d < 1 - n - j. 


Indeed, by the definition of the hypothesis implies that for any I G N, 

(3 G and compact subset K of 


sup 

y'eK,x„>o,i'eR"-^ 


Idlj'^^kiykxr^Ol 

(1 + |^/|)d+i+|/3| 


< oo. 


Thus by the Riemann-Lebesgue lemma, dl^d^,k{y',Xn,x') is continuous and uni¬ 
formly bounded as y' G K, x G R” if d + I + \/3\ < 1 — n, proving the claim. 

Next, recalling the polyhomogeneous expansion 

Af-l 

(37) k = kd-i-j + k^^\ k^^^ G , N = 0,1,2,..., 

j=o 


of k, we claim that 

(38) kd-i-j{y',x) = Ui-n-d+j{y',x) + Qi-n-d+j{y',x) log \x\ + Sj{y',x), 

with Ui-n-d+j £ (^“(R”"^ X R!J:) homogeneous in x and Qi-n-d+j £ (^“(R”^"^ x 
R"') a homogeneous polynomial in x of the degree indicated by the subscript, and 
Sj G C'“(R”“^ X R!(). Combining this with (37) and (36) leads, exactly as in the 
proof of Corollary 5, to the formulas in the theorem for = R!J:, while passage to 
local charts as in the proof of Corollary 7 finally yields the result in general. 

From now on, we thus assume (replacing d — 1 — j byd—1 for brevity) that 
k = kd-i satisfies the homogeneity condition 

(39) k{y'= k^k{y',Xn,i') for A > 1, |C'| > 1. 

Since y' enters — from the point of view of our inverse Fourier transform in — 
only as a parameter on which everything depends smoothly, we will drop it from 
the notation for the rest of the proof. Define l{xn,C') := \C'\‘^k{\^'\xn,C'/\C'\)^ or, 
in the polar coordinates ^' = rQ (r > 0, C G S”“^) on R"""^ 

l{xn,rC) ■■= r'^k{rXn,C)- 

Then I G C'°“(^ x (R^^-^ \ {0})), 


l{xu,0 = k{xn,i’) if l^'l > 1, 

(40) = Vt >0,CV0,x„ >0. 

For brevity, we also denote '■= ^(C^Oi introducing the dilation operator 


hf{x) := f{tx), t > 0, 
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we thus have 

Since k G C'°“(R!Jl) belongs to 5+ as a function of Xn for each fixed , it follows 
that Oil ^ <S+ for any patch function 0, while 


(41) 


h{rC) = 


N-l 

E' 

j=o 


d+j r 


j! 


d^k{Xn,0 


dxi 


+ 0(r^+^) as r \ 0. 


J Xn=0 


Thus in the notation of §4.1, li G Cd(R"'“^). 

Assume first that d > 1 — n, so that li — as well as any It, t > 0 — is integrable 
as the origin, and thus defines a distribution on all of R”“^. For its inverse Fourier 
transform, we have by Proposition 6(ii) 

h G + evi-n-d\og = 

and 

= (x-^S^Ji^ix) = x\-'^-%( — '), Xn > 0. 

\ Xyi ' 

The right-hand side of the last formula is a smooth function on Rl( homogeneous of 
degree 1 — n — d, and from li G it follows that it actually extends smoothly 

to R” \ {0}: namely, if has the polyhomogeneous expansion 


J=0 

then 

OO 

i=o 

We now claim that 

(42) k{xn,x') — lx^{x') =: so{x) satisfies sq G C'°“(R!j:). 

Setting Ui-n-d{xn,x') = lx„{x') we then get (38) (with Qi^n-d = 0), thus settling 
the case of d > 1 — n. 

It remains to show (42). Fix a cutoff function y with x{^') = 1 for |,^'| < 1 and 
/t(C0 = 0 fo^ — 2, and let d = 1 — X be the corresponding patch function. Then 
6k = 91, so 

So = {xk-xiy ■ 

Differentiating under the integral sign in 


{xk^)'^{xu,x') = [ e“ Xi^')kixn,^')d(' 

4|TI<2 

and noting that k G C'°“(R_|_ x R"""^) and k{-,^') G 5+ for each fixed hence 
dZk{Xn,e) ^ X K) for any compact K C R” it is immediate that 

ixky G C°°(R;). 
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On the other hand, from l{xn,rC,) = r‘^k{rXn,C) and the fact (again) that k G 
C'°°(R_|_ X and k{-,C) £ for each fixed C £ 8 ""“^, we see that 

dZl{xn,rO = r‘^^^{dZk)irXn,C) 

with dZ^k bounded on R+ X S” Differentiating under the integral sign in 

{Xinxn,x')= [ e^^'<'x{Z)l{Xn,Zn', 

it again transpires that 

This completes the proof of (42) and, consequently, (38) for the case of d > 1 — n. 

For d < 1 — n, one can reduce to the previous case upon differentiating m times 
with respect to where m = [2 — n — d]: the kernel dZ^k belongs to and 

satisfies (39) with d replaced by d + m > 1 — n. By what we have already proved, 
we thus have 

(43) dZ„k = {dZjr = u^-u-d-m + s,, 

with U G C'°°(R!|:) homogeneous of the indicated degree, and sq G (^“(R”). 

Observe now that if / G C'°“(R!J: \ {0}) is homogeneous of degree s < —1, 
then f{x) dxn (exists and) is homogeneous of degree s + 1 ; if / is homogeneous 
of degree s > — 1, then Z f(x) dxn is homogeneous of degree s + 1 ; and if / is 
homogeneous of degree — 1 , then expressing the difference f{x) — ax^xZ"^, where 
a = /( 0 , 1 ), on the hyperplane = 1 as 

n — 1 

1 ) “ 1 ? /|2 = ^(^0 + \x'?9{x') 

1 ~r \^ \ . ^ 

J = 1 

with fj,g G C'°“(R"'“^) and x a cutoff function on R""”^ (this is achieved simply 
by taking any fj such that Xjfj{x') has the same Taylor expansion around 

x' = 0 as the left-hand side), we get 

n — 1 

(44) /(x) = ^ + XjFj{x) + |x'pG(x), 


where the functions 


Fjix) := 


-X 



G{x) := -^g 
xi 



are smooth on R!J: \ {0} (this is clear for Fj, and follows from (44) for G) and 
homogeneous there of degrees —2 and —3, respectively, so the function 


n — 1 


F(x) := a log |x|Xj / Fj{x) dxn + \x'\^ / G{x)dxr, 


i=i 



22 


M. ENGLIS 


belongs to \ {0}), = / and F{x) — alog |x| is homogeneous of degree 0. 

Integrating (43) with respect to Xn in the above way m times, it transpires that 
k has the form 

m — 1 

(45) Ui-n-d{x) + Ql-n-d{x) log \x\ + So{x) + ^ xivj{x'), 

j=0 

where Ui-n-d £ (^““(RIJl) is homogeneous of the indicated degree, Qi-n-d is 
a polynomial of the indicated degree (in particular, (5 = 0ifl — n — d^N), 
So := (fg dxn)"^so G and Vj are some distributions on R"'“^. Since k 

is clearly a Schwartz function for each fixed Xn > 0, so must be (45), implying 
that in fact Vj G C'°°(R’^“^). Setting Sm ■= Sq + Y^^=o ^ (^““(R”), we thus 

get (38). This settles the case d < 1 — n, and thus completes the proof of the whole 
theorem. □ 

Remark 9. Note that the last proof gives actually a somewhat more precise infor¬ 
mation than (38): namely, if kd-i G R"""^, 5+) is of the form 

kd-i = x^kd+p-i with some p G N and G R”“^, 5+), 


then 


kd-i{y', x) = x'^Ui-n-d-piy', x) + x^Qi-n-d-p(y', x) log |x| + x^s(y', x) 

with Ui-n-d-p G C'°°(R"'“^ X RIJ:) homogeneous in x and Qi-n-d-p G C'°°(R"'“^ x 
R”) a homogeneous polynomial in x of the indicated degree, and s G C'°“(R”“^ x 
R!J:). That is, loosely speaking, if it is possible to pull out a factor of from 
k, then the same is true for k. This is immediate from the obvious fact that the 
inverse Fourier transform with respect to commutes with multiplication by x^. 
It will prove useful in Section 5. □ 

4.3 Singular Green operators. The Schwartz kernel of trace operators of class 
zero can be obtained simply by taking the adjoint of potential operators; the case 
of general class r > 0 then follows easily from (18). It therefore remains to deal 
with singular Green operators. 

Denote by S!( = S"" n the closed upper hemisphere in R"'+^. 

Theorem 10. Let A R”) and denote by y = {y',—yn) the image of y = 

{y',yn) G R” under the reflection with respect to the hyperplane Pn = 0. Then the 
Schwartz kernel k^ of A satisfies 

kA{x,y) = \z\~'^~‘^f(^x',\z\,-^'^ -bG'(x', 2 ;)log|z| 
if d gTj, d> —n; 

kA{x,y) = f[x', l^l, 1^^ G(x',2;)log jz], 

with G vanishing to order —n — d at z = 0, if d G Z, d < —n; and 
kAix, y) = \z\~'^~‘^f(^x', | 2 :|, G{x', z) 
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if d ^ Z, for x,y G R+, where 

z = z{x, y) := {x' - y', Xn,yn) e 

and F G C'°“(R’^-i x R^ x G G x R’"+i). 

In particular, kA extends smoothly up to the boundary of R!J: x R” away from 
the boundary diagonal x = y. 

Notice that an equivalent statement is again obtained upon replacing F{x', \z\, -^) 
by F{x',y', \z\, j^), or even F{x,y, \z\, jly), and similarly for G. Also, by Seeley’s 
extension theorem, one may replace SIJ: by S” or even R"'^^. 

Proof. The idea of the proof closely parallels that for Theorem 8, so we will be 
more brief. By (20), where we again take the kernel g in the y'-ioim. rather than 
the x'-form, 

kA{x,y) = g{y',Xn,yn,x' - y') 

where g denotes the inverse Fourier transform of g{y',Xn,yn,C') with respect to f,'. 
An argument involving the Riemann-Lebesgue lemma again shows that g G x 

R!J:) if d < — n —j, and as before it follows that it is enough to deal with each term 
of the poly homogeneous expansion of g separately. We thus assume from now on 
that g G S'^[~^(R"“^,R"'“^,5_|__|_) satisfies 

(46) g{y',^,^,t^') = t‘^^^g{y',Xn,yn,C') for t > 1, |C'| > 1. 

Omitting again y' (which plays only the role of a smooth parameter throughout) 
from the notation, set for = r(, r G R+, C £ 

l{xn,yn,rC) := r‘^+^g{rXn,rynX)■ 

Tlien I G C'°“(^ X R7 X (R*^-! \ {0})), 

li,Xn, ym ^ ) dXnj Vni C ) if I ^ 1; 

and 

(47) l{^, ^,te) = t‘^^H{xn,ynX') Vt >0X'X 0. 

We further denote 

IrAO :=l{Rcos^,Rsm^,^'), i?G^,i?G [0, f],^' G R^“^ \ {0}; 
in terms of the dilation operator Stf{x) := f{tx), we thus have 

Ir,^ = R ^ 

Since g G C'°°(R+ X R+ X R" ^) belongs to as a function of Xn,yn for each 
fixed which implies that g{Rcos'd, Rsinil,^) belongs to 5_(_ as a function of R 
for each fixed X and D, it follows that Oli^u £ 5+ for any patch function 6, while 

_ -\-k 

^h^rC) = E Pkl (cosd)^'(sind)^[a^^a,^^ff](0,0,C) +0(r^+"+^) 

j+fc<A-l 
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as r \ 0. Thus in the notation of §4.1, G 

Assume first that d > —n, so that — as well as any R > 0 — is inte- 
grable at the origin, and thus defines a distribution on all of For its inverse 

Fourier transform, we have by Proposition 6 (ii) 

h,^ G = 5-"-^ + ev.n-d log = 

(since —n — d < 0), and 

IrA^') = iR-''-^SRhA"'{x') = R—R > 0. 

Since depends smoothly on d, the right-hand side of the last formula is a smooth 
function of {xn,yn,x') = {Rcosd, Rsmd,x') G R+ x R+ X R” homogeneous 
of degree —n — d, and from G it follows that it actually extends 

smoothly to R_|- x R_|_ x namely, if has the polyhomogeneous expansion 

OO / 

iiAx') ~ 

j=0 ' ' 

then 

OO f 

IrA^') ~ 

i=o I* I 

so IrA^') Is C°° in {R,d,x') G R+ x [0, x R"-"^. Arguing as in the proof of 
Theorem 8 , one finally checks that for any patch function 6 on R”“^ vanishing on 
|.^'| < 1 , and X = 1 ~ ^ the corresponding cutoff function, both {xoY {xA 
and, hence, also {xg — xA = {k — A, are smooth on R_|_ x R_|_ x R”“^. Thus, 
setting U-n-diUn, x) = l{xn,yn,x') and restoring the variable y', we have proved 
that for d > —n and g G satisfying (46), 

g{y', Xrr,yn,x') = U-n-d{y', z) + s{y', z), 

z = ixn,yn, x') G R7 X R^ X R"”! =: RlW y' G R^-\ 

with U-n-d £ C'°“(R"“^ X (R+i*).^ \{0})) homogeneous of the indicated degree in z, 
and s G C°°(R^-i x R!(+^). 

Finally, the case of d < —n is again handled by differentiating g with respect to 
Xn (or yn) m times, where m = [1 — n — d], and then integrating m times exactly 
as in the proof of Theorem 8 . The conclusion is that, for any d G R, 

g{y', Xn,yn,x') = U-n-d{y', z) Q-n-diy'^ z) log |z| -b s{y\z), 

(48) -^ 

z = (x^,y^,x')eRA.\y' 

where U G C'°“(R"'“^ x (R!(i*)i^\{0})) is homogeneous in z and Q-n-d G C'°“(R"'“^ x 

Rn+l 

) is a homogeneous polynomial in z of the indicated degrees, and s G (^“(R” ^ x 
R++^)- Putting everything together, the assertion follows in the same way as in 
the proof of Corollary 5. □ 
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Remark 11. Just as in Remark 9, the last proof gives a bit more precise information; 
namely, if a factor of x^y^, with some p, g G N, can be pulled out from g G 
in the sense that 


9 = xlvlG 


with G G 


then it can also be pulled out from the U, Q and s in (48). □ 

Note that in the formulas in Theorem 10, one can replace \z\ and F{x', \z\, |^) G 

C°°(R-i,R 7 , S^) by |x - y\ and F{x, \x - y\, 3) G C°°(R^ x R^ x S^), 
respectively, since the ratio 


(49) 


X - y\^ ^ \x' - + {Xn +ynf 

kP \x'— y'\^ + x"^ + y"^ 


extends to a smooth positive function of (x, \x—y\, |^_|| ) G R” xR+ xS”. Similarly, 

l^l can be replaced by dA{x,y), the distance of (x,y) G R” x R” to the boundary 
diagonal {(x,x) : Xn = 0}. 

The last theorem is again easily transferred from R” to manifolds by means of 
local charts. Recall that for fl a bounded domain in R"" with smooth boundary, 
for e > 0 small enough the mapping 


TT : dQ X (— e, +e) — )• R"", 7r{C, t) = C + 

where is the unit inward normal vector at () G dkl, is a diffeomorphism. Denote 
its image by Ve, and define the map y y: Ve —)■ by 


7r(CT) = 7r(C,-t) 

(the “reflection with respect to 5D”). 

Corollary 12. Let kl be a bounded domain in R"" with smooth boundary, e he as 
above, and A G 0q(D). Then the Schwartz kernel kA of A is (7°° outside Ve /2 x Ve/ 2 , 
while on Ve x Ve it satisfies 

kA{x,y) = \w\~'^~‘^F(x,y, |u;|, -^) + G(x,y)log |r(;| 

V |r(;| / 

if d ^ Z, d > —n; 


kA{x,y) = F(^x,y,\w\,-^'^ + G{x,y)\og\w\, 

with G vanishing to order —n — d at x = y, if d ^ Z, d < —n; and 
kA{x,y) = \w\~'^~'^F(^x,y,\w\,-^'^ +G{x,y) 


if d ^ Z, where 

and F G C'°“(n x D x R^ x 


w = X — y 

G eG°°(nxn). 
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Again, F{x,y, |u>|, |^) could be replaced by F{x, |u;|, |f|) or F(y, |u;|, |^), or ex¬ 
tended smoothly to all of 

Proof. Note first of all that all this holds in the local chart R”, by the last theorem 
(and the remarks after it), since z = (x', y„) is a smooth function of x,y, while 

\z\ is a smooth function of x,y, |rc| and (cf. (49)). The general case now follows 
as in the proof of Corollary 7. □ 

Finally, the Schwarz kernel of an arbitrary singular Green operator of class r > 0 
can be written down easily using the last corollary, (19) and Theorem 8; we omit 
the details. 


5. The Poisson kernel 

Recall that the kernel k{x',Xn,f') G of a potential oper¬ 

ator (17) has the poly homogeneous expansion 

OO 

(50) k{x',Xn,C) ~ '^kd-jix',Xn,f') 

1=0 

where kd-j G ^ ^(R”“^,R"^“^,5_|_) possess the homogeneity property 

(51) kd-j{x', ^,\0 = X^-^kd-j{x',Xn,i') 

for A > 1 and |^^| > 1. It is sometimes convenient — and we will do that from 
now on — to redefine the kd-j by homogeneity on |,f'| < 1 so that (51) actually 
holds for all A > 0 and f,' ^ 0; of course, kd-j will then in general have a singularity 
as = 0 (i.e. will belong only to (^“(R""^ x R_|_ x (R*^"^ \ {0}))), and also the 
expansion (50) will then hold only for |.^'| > 1, i.e. 


(52) 


sup 




k - kd-j 


|^|d— A^+Z— fc—| ck| 


< OC 


for all k,l, N > 0, multiindices a, (5 and compact subsets K of R""”^. 

We have seen in the proof of Theorem 8 in §4.2 above that each homogeneous 
component 


(53) Ui-n-d+j{x ,z) + Qi-n-d+j{x ,z)\og\z\ 

in the Schwartz kernel kxix, y') = k{x', z), z = (x„, x' — y'), of the potential oper¬ 
ator in question (with Ui-n-d+j, Qi-n-d+j homogeneous of the indicated degree 
in z — the latter a polynomial — and smoothly depending on x') actually arises, 
modulo an error term smooth on all of R”“^ x R”, precisely as the inverse Fourier 
transform with respect to of 


kd-j{x',Xn,(')0{f,') 

with some patch function 6. We use this correspondence to obtain some extra 
information about the functions Ui-n-d+j,Qi-n-d+j in (53) in the particular case 
of our Poisson operator K. 
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So let us consider our Dirichlet problem 

A/= 0, f\9Q = u 

on a bounded domain C R” with smooth boundary dQ. We want to compute 
the components (53) of the corresponding Poisson kernel K{x, C) for x, C near some 
boundary point a G dQ. Without loss of generality, we may assume that a is the 
origin and that the tangent hyperplane to at a is given by = 0. In a small 
neighbourhood of the origin, 0, and will then be given by 

(54) Q = {x : Xn > (/>(x^)}, = {x : Xn = 4>{x')}, 

for some smooth function cj) on satisfying (/)(0) = 0, Vi;i)(0) = 0. In the same 

neighbourhood, we can therefore use the simple map 

4'(x) = {x’,Xn - (t){x’)) 

as a local chart mapping (a piece of) onto (a piece of) R” . The Laplace operator 
A on n transforms under T into the second order operator 

£ = (I + |V(/>P)9nn - {^4>)dn “ 2 ^ (9fc(/>)9fcn + ^ 

k^n k^n 

= (1 + |V</>|2)a„„ - {Acj))dn - 2V(t> • Vdr, + A' 

on R!j:, where we write for brevity dk = g§- and djk = djdk, and A', V' denote the 
Laplacian and the gradient with respect to x' (the prime being omitted if there is 
no danger of confusion). 

Let us now look at the potential operator K which solves the transformed prob¬ 
lem, 

Cf = 0 on R”, /(-jO) = u on R”“^. 

We know from Section 3 that K is of the form (17) with k G (R"""^, R”“^, 5+) 
where d = 0. The boundary condition clearly forces (cf. (24)) k{x',0,^') = 1 for all 
x' and , or, in terms of the homogeneous components (50), 

(55) ko{x',0,^') = 1, k-j{x',0,^') = 0 foT j > 0, Vx^,G R””^. 

On the other hand, CK = 0 means that 

(56) £[6*""/c(x',Xn,C0] = 0 Vx',x„,^' 

(£ is being applied to the x',x„ variables). Now 

e-i-<'c[e^-<'k{x', xn,e)] = (1 + I - (A0)a„fe - 2V</) • {iCOnk + Vd^k) 

+ {A'k + 2iC ■Vk-\C'\‘^k). 

Substituting for k the polyhomogeneous expansion k-j and comparing terms 
with the same degrees of homogeneity, (56) yields the system of equations 

(57) Mk-j = TZk-j+i + Zk-j+2, j = 0,1,2,..., ki = k2.= 0, 
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where Ai,TZ,Z are ordinary differential operators in Xn- 

Mg = -(1 + \V(l)f)dr,ng + m' • ^4>)dng + 

TZg = -{A(j))dng -2V(j)- V’dng + 2i(^' • V’g), 

Zg = A'g. 

Also, k-j {x', Xn,i') need to decay rapidly as Xn —)• +oo for each fixed x', (because 
they belong to 5_|_ with respect to Xn)- From this it is possible to solve (57), with 
the boundary conditions (55), recursively with a unique solution at each step and 
obtain all the k-j. 

For j = 0, (57) just reads Mko = 0, so ko is a linear combination of '> 

and \ where 

^ ,,, _ • v <))± + lv<(>P) - le • v </.|2 

l + | V ^ 6|2 

The requirement of rapid decay in Xn means that only enters, and (55) then 
yields 

ko{x',Xn,0 = 


Next, for j = 1, (57) reads 


Mk-i = TZko- 


Note that r/_|_ is homogeneous in of degree 1. A moment’s computation therefore 
reveals that the right-hand side of (58) has the form 

'kiko = Fi{x',C)ko + XnF2{x',C)ko 

where Fi,F 2 are homogeneous in of degree 1 and 2, respectively. Solving (58), 
one finds that k-i has to be of the form G{x' ,^')ko+Go{x' ,^')xnko+Gi{x' ,^')x‘^ko. 
The boundary condition (55) forces G = 0, and thus 

fe_i(x',x„,eO = [Goix',e)xn + G^ix',e)xl]F^-^+^^'’^'\ 

with Go,Gi homogeneous in of degree 0 and 1, respectively. Continuing in 
this way, it transpires that for all j > 0, 


k.,{x',x^,e) = 


for some Gqj{x' homogeneous in of degree q — j- 

The contribution from k-j to the Schwartz kernel kxix,!/) of K (i.e. to our 
Poisson kernel K{x,Q in the local chart T) is thus given by the inverse Fourier 
transform with respect to of 
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evaluated at x' — y'. (For j = 0, this is to be interpreted just as \) 

In particular, at the point x' = 0 that we are interested in, we have V(/>(0) = 0 so 
and we need the inverse Fourier transform with respect to of 

2j 2j 

(60) =: 

q=l q=l 

(we omit the argument x' = 0 in Gqj for brevity). Analogously as in §4.2, one shows 
that for g — j > 1 — n, i.e. when F)jj,a;„(^0 = I is integrable at the 

origin, Fqj^x^{x') = Fqj,i{^) is homogeneous on R” of degree 1—n—g+j; 

secondly, as F’gj,! C Cg-jL whence C g^-n-q+j^ Fqj^Xr,{x') is smooth 

in {x',Xn) G R” and extends smoothly to R” \ {0}; thirdly, that inserting the 
patch function 9 introduces only an error term smooth on all of R”; and finally, 
differentiating and then integrating back m times with respect to where m = 
[2 —n —g+j], that all the facts just mentioned remain in force also for q — j < 1 — n. 
Thus the inverse Fourier transform of (60) with respect to has the form 


2j 

( 61 ) -q+j{x) + Ql-u-q+j{x) log |x|] + s(x) 

q=l 

with Ui-n-q+j , Qi-n-q+j smooth and homogeneous of the indicated degrees on R 4 :\ 
{0}, with Qi-n-q+j a polynomial, and s G C'°°(R”). Again, for j = 0, the sum in 
(61) has to be replaced just by the inverse Fourier transform with respect to of 
g-ainIS 1^ which is well known to be equal to (see e.g. [22], p. 247) 


(62) 


CnXn 

(|x'P + x2)"-/2 ’ 


— 


r(f) 

yj-n/2 


For j > 0, the log-term in (61) appears only for 1 — n — g -|- j > 0; since q > I, 
it therefore occurs only for j > n. Summing over j, it thus transpires that the 
Poisson kernel K{xX) = k{x',Xn,x' — y') in local chart satisfies 


k{y\xn,x') = j^F{y',\x\,^\ -L XnG(y', x) log |x| + s{y\x), 


with Cn as in (62), F G C'°°(R"“^ x R+ x S”“^) and G,s € C'“(R”“^ x R”). 
Since, from its very definition, K{x,y') vanishes for x^ = 0 (and x' / y'), it must 
be possible to pull out a factor of x^ also from s{y',x), and then merge the latter 
into the first summand. Altogether, we have thus arrived at the following more 
precise version of Theorem 8 for the particular case of the Poisson kernel, which 
we have stated as Theorem 1 in the Introduction. 

Theorem 13. (= Theorem 1) Let LI be a bounded domain in R"" with smooth 
boundary. Then the Poisson kernel K{x,C) for the Dirichlet problem (2) can he 
written in the form 




X - Cl, 


a: - C \ 

X - Ck 


-LG(x,C)|a: 


Criog|x 


Cl 


i^(x,C) 
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where d G d > 0 on and d{x) = dist(x,9i7) for x near dQ, Cn is as 

in (62), F G X R7 X F(-,0,-) = l, and G ^ (U x dO). 

We conclude this section by showing that the log term is indeed present in 
general, although the explicit formulas 

K(x C) = 

^ 2 \x- Ch 

for n = {x G R” : |x| < 1}, and 

IX( ^nXn 

A (x, C) = 1 - TT- 

for n = R”, show that it is absent for the ball and the half-space. To that end, 
consider the special domain 

f] = {y G R" : 2/„ > <h(|y?)} 

for some smooth function <f> on R+; i.e. the special case of (54) where (p{y') = 
‘h(ly'p) depends only on |y'|. The computations above then simplify somewhat; 
for instance, one has the explicit inverse to the operator A4, namely, the unique 
solution to 

Mg = x(^e“"’'+ 

(m = 0,1, 2,...), rapidly decreasing as x„ —)• -|-oo, is 


9 = 


ml 

(_y)m+2 


m + 1 

E 


i=i 


{ — XnVy e*a;„77+ 

j\im+2-j 1 + |V0|2’ 


where v '=9+ — V-- Using all this, one arrives at the following formulas (checked 
on a computer) for n = 3 and x' = 0: 


ko{Q,Xn,C) 

k-i{Q,Xn,C) 

A:_2(0,Xn,CO 

/c_3(0,Xn,^') 


- g 1^ I 


= x„(l-x^|C'|)$'(0)e 

= (^ + - 3|{'k + kM) 


2|e|3 

19$'(0)3 - $"(0) 


2|el 

5$'(0)" 




4>'(0)= 


^"( 0 )) 


and, hence, after a small computation^ 


F((0,Xn),x') 


( 63 ) 


t t^(2-3t^)4>^(0) 

27r|xP 27r|x| 

t(l - f2)(l + 12^ - 15t4)4>'(0)2 

A$"(0) - 24>'(0)3) 

H--- -\x\ log |x| -|- (higher order terms). 


^See (70) below for a hint. 
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where we have set for brevity t = Xn/\x\. (All these formulas relate to the local 
chart by obtained from Q via the map {y',yn) {y'ly-n — = {x',Xn).) 

Thus the log term is present whenever d>"(0) — 2<h'(0)^ 7 ^ 0. Note that the last 
expression vanishes if Q is the half-space (i.e. = 0), or if n is the ball of radius 

R tangent to ?/„ = 0 at the origin (i.e. ^{t) = R — VR^ — t, so that <h'(0) = 
<h^^( 0 ) = in complete accordance with the fact that there is no log term in 

K{x,C) in these cases. 

Finally, we remark that the log-term is always absent in dimension n = 2. 
Indeed, since harmonic functions in C = are preserved by composition with 
holomorphic maps (this is no longer true in C*” = R^*” when m > 1 , nor has any 
good analogue in odd dimensions), while any simply connected domain in C with 
smooth boundary can be biholomorphically mapped onto the disc (by the Riemann 
mapping theorem), and the Poisson kernel of the disc has no log-term, the Poisson 
kernel of the original domain has no log-term either. (The argument extends in 
fact also to multiply connected domains with smooth boundary, in view of the local 
character of the boundary TDOs.) Thus the examples above for re = 3 are probably 
the simplest domains for which the log-term is present. 

6. The harmonic Bergman kernel 

The singularity of the Schwartz kernel of a singular Green operator G can again 
be obtained from the asymptotic expansion of its symbol 

00 

(64) 5 ~ X] 

j=o 

into the homogeneous components satisfying 

(65) gd-j{x',^,i^,X^') = X'^~^^^gd-jix',Xn,yn,^') 

for A > 1, |.^^| > 1; as with potential operators, we will again assume from now on 
that (65) actually holds for all A > 0 and / 0, at the expense of having gd-j 
singular at = 0 and the asymptotic expansion (64) holding only for > 1, 
analogously as in (52). The contribution to the Schwartz kernel 

( 66 ) kG{x,y) = g{x',Xn,yn,z'), z'= x'- y', 

of G from gd-j is again given by the inverse Fourier transform with respect to of 

9d—j(x : Xn: yn} ^ )^(? ) 

for some patch function 9. To obtain the gd-j more explicitly, we use the result for 
the Poisson operator from Section 5 together with the composition rules reviewed 
in §2.3. 

Recall that our singular Green operator of interest, the harmonic Bergman pro¬ 
jection, is given by (31) 


G = KA“^K*, A:=K*K, 

where K is the Poisson extension operator from Section 5. By the general rules 
for the Boutet de Monvel calculus, A is a TDO on the boundary of order —1, and 
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G is a singular Green operator of order 0 and class zero. Let us now compute the 
symbol of G. In terms of the symbol k{x',Xn,^') of K (in a local chart by R”), 
the adjoint K* is a trace operator of class zero with symbol k*{x',Xn,^') given by 
(cf. (25)) 

k*{x',Xn,C) ~ 

a\ ^ 

a 

This formula is well adapted with respect to the homogeneous expansion (64): 
namely, one gets (remember that K is of degree d = 0) /c Ejlo with 

ra-\-\oi\=j 


We have seen that k-j are of the form (59) (where for j = 0 the sum YllLi Gqjx'l^ 
has to be interpreted as 1). The last formula and a quick check show that k*_j are 

of the same form, only with replaced by ; in particular, 

k* = 

Next, A = K*K is a TOO on dVL (= in the local chart) with symbol given 

by (22) 

s(a:',^') = X] ^ J D^'k*ix',Xn,Od^'Hx',Xn,^')dXn. 


This is again well-behaved with respect to the polyhomogeneous grading, and using 
the fact that 





k\ 

[if]- — Z7?+)^ + ^ ’ 


we get s ~ Xljlo with 


1 

s-m-i= ^ D‘^,k*_jd^,k-idxn, 

l+j + \a\=m ' ^ 

and, in particular, s_i = lj{ir]- — i??+). The symbol p of the inverse A“^ therefore 
has the asymptotic expansion p ~ Yl^=oPi-ji given recursively by 

Pi = —, Pi-m = - —DZpi_jd'^,S-i-k for m > 0, 

s_i s_i a! ^ 

k-\-j-\-\oL\=m,j<m 


by (15); in particular, pi = i(ry_ — 1 ?+). 

By (22) again, we get that KA“^ is a potential operator of degree 1 with symbol 
X ~ Xi-j with 


vi-^ = ^D'^,k-jd^,pi-i. 

j + \a\+l=m 

For the terms with I < m, (59) and a routine check imply that they are again of 
the form (59) (with m in the place of j) multiplied by a function homogeneous 
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in of degree 1 — Z; while for I = m, the corresponding term is just 
Altogether, we thus see that vi-m has the form 

2m 

(67) ui_^ = 

q=0 

with Fmq homogeneous in of degree 1 — m + q and Fmo = Pi-m- (Note that 

(67) is similar to (59), except that the sum starts from q = 0.) In particular, 
Vi = i{ri- — r/+)e*'^"^+. 

Finally, using ( 22 ) one more time, we get that KA“^K* = G has symbol g ~ 
J2T=o9i-j, with 

9l-m= ^D'^,Vi-j{x',Xn,C)dx'k-lix',yn,^')- 

j + l + \a\=m 

In particular, gi = i{g- — , while for general m we again have from 

(59) (for k*) and (67) 

2m 

(68) 51-rn = Fmrqix',eWne^^^-^+-y-^-\ 

r,q=0 

with Fmrq homogeueous of degree 1 — m + r + g in and F^oo = Pi-m- 

The contribution to the singularity of the Schwartz kernel ( 66 ) of G — i.e. to the 
harmonic Bergman kernel H — is given by the inverse Fourier transform with 
respect to evaluated at x' — y', of ( 68 ) multiplied by a patch function 
Going through the proof in §4.3 (analogously to the argument after (60) in the 
preceding section), it transpires that the latter equals, modulo smooth functions 
on X 


(69) 




U —Yi-\-rn — r — q{p^ ' ,z) + Q 

— n-\-m — r — q{p^ ',z) log\z\ 


Z = {x! - y',Xn,yn) G 


for some U-n+m-r-qix' , z), Q-n+m-r-q{x' , z) smooth ou R*^"^ X (R”"*"^ \ {0}) and 
homogeneous in z of the indicated degree, with Q a polynomial in z. 

In particular, at x' = 0 , we have ig- = —ig+ = |^^|, so gi = I, 

and 


(70) 


9i = 


2— 

V 

_ Cnt 

L dt \ 


t=x„+y„ dt {F + X' 2)"-/2 


t—Xfi ~\~yn 


with Cn as in (62), which shows that the leading term of kc{x, y) in the local chart is 

^ {n-l){xn+ynY-\x'-y'\^ 

''"[(Xn+yn)2 + |F-y'|2]-/2+l- 


In terms of the boundary distance d{x) and the reflected point y = (y', —yn), this 
becomes 

^ n{d{x) + d{y))^ - \x - y\‘^ 

\x-y\n+2 

Since the differential of our local chart T from Section 5 equals the identity at the 
origin, the last expression remains in force also back in fl. We have thus arrived 
at the following slight strengthening of Corollary 12 for the particular case of the 
singular Green operator G = KA“^K* of order 0 and class zero, which has been 
stated as Theorem 2 in the Introduction. 
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Theorem 14. (= Theorem 2) Let LI be a bounded domain in R"" with smooth 
boundary. The harmonic Bergman kernel H{x,y) of id is C°° for {x,y) away from 
the boundary diagonal, while near diagSO it can be written in the form 

(71) H{x,y) = ^J‘y^^^ F[x,y,\x -y\, |^_ ~| ) + G{x,y)\og\x - y\, 


where G G F G C°°(OxOxR+x F{x,x,t),v) = n{v,Vd{x))‘^-l 

for X G dil, and Cn is as in (62). 

Note that, in particular, the leading term of the restriction i7|aoxao of H to the 
boundary is simply —2cn|x — y\~'^. 

Remark. In order to obtain an expression for H{x, y) like (71), valid on all of H x fl, 
not only near the boundary, one possibility is as follows: 


(72) 


H{x,y) 


2cjj 

v{x,y)\^ 


F 



v{x,y) \ 
v{x,y)\J 


+ G{x,y)log\v{x,y)\ , 


where G G C'°“(0 xLl), F £ C°°{Ll x x R+ x S’"+^), and v: Q x Q ^ R’"+i is 
any smooth function such that 


v{x, y) = {p{x) - p{y), d{x),d{y)) 

for {x,y) in some small neighbourhood of diagSfl, where p{x) is the point of dil 
closest to X (p is well defined in a small neighbourhood of dil, and such v exists by 
Seeley’s extension theorem). One can further replace \v{x,y)\ in (72) by d/\{x,y), 
with any d/\ G G°^{il x 17) such that d/\ > 0 on 17 x 17 \ diag517 and d^{x,y) = 
dist((x, y), diag 917) for {x,y) near diag917, etc.; cf. the remark after the proof 
of Theorem 10. □ 

We conclude this section again by showing that the log term is present, i.e. G ^ 0, 
in general, although this is not the case for the familiar situations like the ball 
17 = {x G R” : |x| < 1} with 

Cnin- 4)|x|^|y|^ + (8x • y - 2n - 4)|xp|yp + n 

2 (l-2x-y + |x|2|y|2)-/2+i 

or the upper half-space R” with 


ut ^ o n{xn + ynf-\x-y\^ 

- 

To this end, we use once more the domains from Section 5 

17 = {y G R" : y„ > <h(|y?)} 


with some G C°“(R+) and n = 3. We again pass to the local chart R!( by means 
of the map 'k ; {y',yn) {y',yn — ‘^’dkl^)) = {x',Xn). We will consider H{x,y) 
only at x' = 0 and Xn = yn = 0 (that is, we exhibit a log-term even among the 
terms (69) with y = r = 0). Now since F^oo = Pi-m, we have that, modulo smooth 



POISSON AND HARMONIC BERGMAN KERNELS 


35 


functions, H (0, {y', 0)) = ^(0,0,0, —y') is the inverse Fourier transform with respect 
to evaluated at —y', of 


5(0,0,0,e') =p(0,O- 

Combining the formnlas for k*, s and p above with those for fco, k-i, k- 2 , k -3 in 
Section 5, we obtain after some calcnlation (again verified on a computer) 


1 <h'(0) , 5<1>'(0)2 , 5<h"(0) , 

? ) ~ TT7T77 - TTTLT + + • • • 


2iei lep 4|^'|3 2ier 


and 


p{o,e) ~2|c'i+44>'(o)- ^ + .... 


1^1 


l^'l 


Consequently, 


^ 3 , p{0, ^ - C3 (24." (0) - 44' (0)3) log 1 9 ' I 

+ (higher order terms). 

Hence the log-term is nonzero as soon as 24>''(0) — 44>'(0)3 ^ 0. Note that, again, 
the last expression vanishes for the upper half-space (<h = 0 ) as well as for the ball 
(4>(t) = R — y/R^ — t, <h'(0) = 4>"(0) = 4 ^), as it should. 

7. Concluding remarks 

7.1 Weighted harmonic Bergman kernels. For a positive weight function 
w on H, one can consider the subspace of harmonic functions in 

the weighted space L‘^{VL,w); if w~^ is locally bounded (in particular, if w is 
continuous), the standard argnment shows that point evaluations are continnons 
on tc), and uniformly so on compact subsets, so that 7y^arm(^)actu¬ 

ally closed in L^(i4, w) and possesses a reproducing kernel — the weighted harmonic 
Bergman kernel H^{x,y). Of course, for not to be identically zero, the space 
-^harm(^>^eeds to be nontrivial; this is the case e.g. whenever w is integrable 
over 0 (since the constant functions then belong to 

Much of what has been said above for the unweighted case remains also true 
in this weighted situation; in particular, the adjoint of the Poisson operator K with 
respect to the weighted inner product on H is just K*rt; (i.e. the unweighted adjoint 
K* preceded by multiplication by w), and the operator 

(74) := K*n;K 
is related to the weighted Bergman kernel by 

(75) H^{x,y) = {A~^K^,Ky)dn, K^ = K{x,-), 
and to the weighted harmonic Bergman projection 11,^ by 


n,, = KA-3K*n;. 
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The Schwartz kernel of however, is not equal to H^,, since the Schwartz kernel 
is taken “with respect to the Lebesgue measure”, while Hyj must be integrated 
against the weight w] for this reason, Hy^ is actually the Schwartz kernel of 

G := = KA-^KT 


In particular, if is bounded with smooth boundary and w is of the form 
(76) w{x) = 


with d G d > 0 on Q and d{x) = dist(x,9fl) for x near a > —1 

(to make w integrable), and g G C°°{Q), then the operator (74) belongs to the 
Boutet de Monvel calculus; namely, is a TDO on of order —a — 1 and with 
nonvanishing principal symbol (see [8]) 

T{a + l)e^G') 

(2|e'|)“+i ■ 


Since K is a potential operator of order 0 while K* is a trace operator of class 
zero and order —1, it follows that G is a singular Green operator of class zero and 
order a, with principal symbol (in the local chart from Section 5, at x' = 0) 


(78) 


9a-\-lix j Xn, yn, ^ ) 


(2|^'|)"+ie-(2:r.+?;r.)|?'| 

r(Q; + l)esG') 


It follows that Hyj{x, y) is given by Theorem 10 (and Corollary 12) with d = a, and 
the leading term (in the local chart by R-lfl) can be obtained as the inverse Fourier 
transform with respect to evaluated at x' — y', of (78). (For a integer, the latter 
can be computed explicitly as 


ga + l{0,Xn,yn,x') 


/ [e *l^'l]^(x') 

V dt) T{a + l)e3G') 


t=X„+yn 


Cn / + ^ t 

F(q; + l)es(^') V dt) (t^ + 




For general a > —1, 5q,_|_i( 0, x^, 2/n, 2:0 is given by the hypergeometric function 


2 °F(^)F( ^+°-i ; 


TT 


,/2r(^)F(c, + l)eSG')(x„+y„) 


n-\-a — l 


2 F 1 


n-\-a n-\-a — l 

2 ’ 2 
n — 1 
2 


J\2 


{Xn + Vn)^ 


as can be seen by expanding into power series and integrating term by term.) 

7.2 Log-terms in dimension 2. We have seen in Section 5 that the Poisson 
kernel of a smoothly bounded domain II C = C has no log-term, due to the fact 
that this is the case for the disc and the Riemann mapping theorem. Computations 
indicate that also the harmonic Bergman kernel H{x,y) has no log-term for such 
domains; or, using again the Riemann mapping theorem, the weighted harmonic 
Bergman kernel Hyj on the disc with weight w of the form w = where / is a 
conformal map, has no log term. It would be interesting to know if this is indeed 
true, and why. (For general weights w on the disc, Hy, has log terms; an example is 
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w{z) = 2—jzp, wheni7^(x,y) = 2Re[(l——(1—21og(l—xy)+F(xy)], 
with F continuous on the closed disc.) 

For our sample domains {x G : X 3 > <h(|x'p)} in Sections 5 and 6 , it was 
also the case that the conditions for the presence of the log-terms in K{x,C) and 
in H{x,y) were the same — namely, <h''(0) — 2d>'(0)^ / 0; suggesting that perhaps 
the log term is present in K{x, Q if and only if it is present in H{x, y). The present 
author has no idea whether this is indeed true, nor does he know whether there 
are other domains than the ball and the half-space for which the log-term in either 
iF(x,C) or H{x,y) would be absent. 

7.3 Harmonic Sobolev-Bergman kernels. Our methods apply also to the re¬ 
producing kernels of Sobolev type spaces of harmonic functions. For instance, 
consider the space 

^harm(^) = harmouic on n : u, Vu G L^(n)} 
with the standard norm 


ll^lli •— Iloilo + ^ ll'9ju||^. 

i=i 

Note that the function djU is harmonic if u is; thus, for / G = 

K.RjU, where 

Rj = -fdjK = A-^K*djK 

is a TDO on dQ of order 1, by Boutet de Monvel’s calculus. Then 

n 

\\Kf\\l = {Kf,Kf)n + J2{KRjf,KR,f)n = {Tf,f)dn, 

i=i 


where 

n n 

(79) T = K*K + R*K*KRj = A + ^ R*ARj 

j=i j=i 

is a positive self-adjoint TDO on d^l of order 1, with principal symbol cr{T) = 
(j(A) \a{Rj)\‘^. One can show that 


i=i 


1 

2cj(A)2’ 


thus cr{T) = icj(A)“^ > 0, so T is in fact elliptic. 

Consider now, quite generally, the Hilbert space of harmonic functions on Q 
obtained as the completion of {K/ : / G C'°“(9H)} with respect to the norm 


K/|||:= (T/,/), 


for some positive self-adjoint elliptic TDO T on dO.. For the corresponding repro¬ 
ducing kernel Hxix^y) = Hx^xiv), we then get 


{f,K^)dn = K/(x) = {Kf,HT,x)T = {Tf,jHT,x)dn 
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by the definition of the T-norm. Hence = TjHt^x: or 

(80) -fHT,x = T-^Kx 
and 

HT{x,y) = K'jHT,x{y) = {jHt^x, Ky)dfi = {T~^Kx, Ky)dfi, 

generalizing the formulas (6) for T = K*K and (75) for T = K*tcK. Finally, Ht is 
the Schwartz kernel of the operator G given by 

GKf{x)= [ Kf{y)HT{x,y)dy 
Jo. 

= {Kf, HT,x)n = (K/, KT-^Kx)n by (80) 

= {T-^K*Kf,K,) 9 n = KT-^K*Kf{x), 

that is, 

(81) G = KT-^K*. 

For the harmonic Sobolev space above, with T given by (79), the Boutet 

de Monvel calculus shows that G is a singular Green operator of class zero and 
order —2, and the corresponding harmonic Sobolev-Bergman kernel Hx is thus 
described by Theorem 10 and Corollary 12 with d = —2. 

For the analogous harmonic Sobolev spaces IFharm(f^) of higher integer order k, 
the corresponding operator T becomes 

k n 

r=E E 

m=0 jl,...,jrn = i 

a positive self-adjoint elliptic TDO on did of order 2k — 1, so G is a singular Green 
operator of order —2k (and class zero) and Theorem 10 and Corollary 12 again 
apply. 

The operator T covers also the situation of weighted harmonic Bergman spaces 
in §7.1 (then T = K*rcK = A^), and one can even combine the two and consider 
weighted harmonic Sobolev spaces with respect to weights w of the form (76). This 
is important due to the fact that Tharm(^> "'^hh w as in (76), is known to actually 
coincide with and similarly coincides with bFh~“^^(n) 

(in both cases, with equivalent norms); varying k and a, one can thus cover the 
whole range of bFhg^i.ni(^) with any real s (without resorting to interpolation, which 
is usually used to define for non-integer s). The corresponding reproducing 
kernels are still susceptible to the treatment as above. 

7.4 Analytic continuation. It has been shown that for weights of the form (76), 
the corresponding holomorphic Bergman kernels Ba{x,y) can in fact be contin¬ 
ued analytically in a to a meromorphic function in the entire complex plane [18]. 
The same idea — in fact, a much simpler variant of it — works also for our harmonic 
kernels here. 

Namely, fix d,gG G“(H), with d > 0 on and d{x) = dist(x,9n) near dil, 
and let 

(82) Waix) := d{x)°‘e^^^\ a > —1. 
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The associated weighted harmonic Bergman kernels are then given by 

the formula (75) 

Hy,ix,y) = {A-^K,,Ky), 

where 

Aq := = K*WaK. 

are positive self-adjoint elliptic d'DOs on di} of order —a — 1. Now, first of all, 
Proposition 10 and Remark 11 in [18] assert that Ka/T{a -|- 1) can be analytically 
continued to all a G C as a holomorphic family of TDOs (which we still denote A„) 
of order —a — 1 (with the principal symbols still given by (77), only with the factor 
r(a -|- 1) omitted). It therefore remains only to deal with the invertibility of Aq,. 
On the other hand, A = K*K is a positive self-adjoint operator, hence we can form 
its complex power A“ for any a G C by the spectral theorem; as A is an elliptic 
TDO of order —1, the familiar theory of Seeley [39] tells us that A“ is actually a 
TDO of order —ol with principal symbol (t(A)“ = (2j^'j)““. It follows that 

F(a) :=r(a + l)“^A-“Ao^A„, a G C, 


is a family of TDOs on cAI, holomorphic in the usual sense (see §2.3 in [18]), 
of order 0, and with principal symbol equal to (cf. (77)) 


1 


r(Q; -|- 1 ) 


(2ie'i)^ 


2]^'] r(a + l)es(^') 


og{x') 


( 2 |el) 


a+l 


= 1 . 


Thus F{a) — / is a holomorphic family of d'DOs of order —1, hence, as dO. is 
compact, of compact operators; furthermore, T(0) = A^C^Aq = /. A theorem of 
Gohberg ([24], Chapter I, Theorem 5.1) implies that there exists a set U of isolated 
points in C such that F{a) is invertible for a , and F{a)~^ is holomorphic on 
C\U with poles with finite-rank residues at points of U. Consequently, 

Y{a) := r(a + l)-iF(a)-iA-“A-i 


is a holomorphic family of TDOs of order a-\-l on C\U, with poles with finite-rank 
residues at the points of U, which coincides with A“^ for a > — 1. Thus 


Hy,{x,y) := {Y{a)K^,Ky) 

gives, for any x,y £ Q, a holomorphic function on C\U, with (at most) poles at the 
points of U, which coincides with the weighted harmonic Bergman kernels 77q,(x, y) 
for a > —1, proving our claim about the existence of the analytic continuation. 

Since the operators Ga '■= KT(q;)K* still belong to the Boutet de Monvel cal¬ 
culus — being singular Green operators of class zero and order a — the boundary 
behaviour of the analytically continued kernels Ha{x, y), a G C\U, is still described 
by Theorem 10 and Corollary 12. 

Finally, let us remark that for Q = {z £ C : \z\ < 1} the unit disc in = C and 
w{z) = w{\z\) a radial weight, it is elementary that the harmonic Bergman kernels 
are, up to the constant term, just the real parts of the corresponding holomorphic 
Bergman kernels: 


Hy,{x,y) = Byy{x,y) -L Byy{y,x) - Byy{0,0). 

Taking in particular w = Wa, this implies that the pole-set U above will be the 
same for the holomorphic and the harmonic kernels. In particular, it follows that 
the pole-set U can assume the various bizarre forms described in §7.2 of [18]. 
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7.5 Logarithmic weights. The results of this paper might extend also to bTDOs 
which are not classical but log-classical (or log-polyhomogeneous), in the sense 
of allowing symbols with the more complicated expansions 

OO ^3 

k ~ ^ ^ ^ ^ , Xn , ^ ) 

j=0 m=0 

where 

kj3n{x', AO = X‘^-H^OgXrkj3n{x',X33,e) 

for l^'l = 1 and A > 2; and similarly for the trace and the singular Green symbols. 
In the language of §7.1, this would lead to the description of e.g. the boundary 
behaviour of weighted harmonic Bergman kernels 77^(a;, y) with weights w of the 
form 

OO 

w{x) ~ (i(x)" EE d{xy (log d(x))™'e®"*^' 

j=Q m=0 

(where Mq = 0 , a > — 1 , gjm G as has been done for the holomorphic case 

in [17]. 

7.6 Berezin transforms. The Berezin transform associated with a weighted holo¬ 
morphic Bergman kernel By 3 {x,y) is the integral operator on 11 defined by 

Byofix) = B33 ,{x,x)~^ [ f{y)\B333{x,y)\‘^w{y)dy. 

Jn 

Note that by the reproducing property, fixes holomorphic and anti-holomorphic 
functions. For weights w = of the form (82), the asymptotic expansion as a —)• 
-|-oo of the Berezin transforms Bw^ = Ba plays crucial role in quantization of Kahler 
manifolds (the Berezin and the Berezin-Toeplitz quantizations, see e.g. [1], [38]). 
It would definitely be of interest to understand the behaviour as a ^ -boo of the 
similarly defined Berezin transforms also in the harmonic, rather than holomorphic, 
setting studied in the present paper. While there are a few results available in some 
special situations, the general case remains unclear. 

7.7 Curvature invariants. In the holomorphic case, the boundary values of 

the functions a,b in Fefferman’s expansion (1) involve interesting biholomorphic 
invariants (see e.g. Graham [25]). For the harmonic case, it is clear from the proofs 
above — and, in fact, is quite immediate from the pseudolocality of bTDOs — that 
the singularity of K{x,() as some boundary point x = ( = a G dQ, or of H{x,y) 
at x = y = a G did, is determined completely by the jet of dil at a. To determine 
how exactly the Taylor coefficients of the F{C,,r,v) G C°°{dil x R+ x S”“^) and 
G{x, C) G X (912) in Theorem 13 at C = a, r = 0 and x = C, = a, respectively, 

or of the F{x,y,r,u) G C'°°(II x II x R+ x S"""^) and G{x,y) G G°°{il, x 17) in 
Theorem 14 at x = y = a, r = 0 and x = y = a, respectively, depend on this jet — 
e.g., whether they depend just on some curvature invariants of dil at a — seems 
to be rather difficult. (A hint is given by the formulas (63), (73) in our examples 
in Sections 5 and 6.) Note that in the holomorphic case, there are many maps that 
preserve holomorphic functions, i.e. functions annihilated by d (all biholomorphic 
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mappings); whereas the only maps preserving harmonic functions, i.e. functions 
annihilated by A, in R"" with n > 2 are just the rigid motions and dilations. 
(See e.g. [3], p. 44.) Consequently, there is no hope of bringing dVt. into some 
simpler “canonical form” while preserving its harmonic functions (like the Chern- 
Moser normal form in the holomorphic case). For this reason, it is quite likely that 
one will need the complete information about the jet of and not just e.g. the 
curvature or similar combined quantities, to determine the singularities of the F 
and G above (or, equivalently, of K{x,Q) and H{x,y)). 
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